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Non-Abelian extensions of fluid dynamics, which can have applications to the quark-gluon plasma, 
are given. These theories are presented in a symplectic/Lagrangian formulation and involve a fluid 
generalization of the Kirillov-Kostant form well known in Lie group theory. In our simplest model 
. the fluid flows with velocity v and in presence of non-Abelian chromoelectric/magnetic E"/B" fields, 

' the fluid feels a Lorentz force of the form QaE" + (v/c) x QaB°, where Qa is a space-time local 

I non-Abelian charge satisfying a fluid Wong equation [{Dt -f v ■ D)Q]a = with gauge covariant 

^S) ■ derivatives. 
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I. INTRODUCTION 



When different species of particles compose a fluid, it may be that the collective variables describing the fluid 
(density p, velocity v in an Eulerian formulation), and the dynamical equations that govern them, can reflect the 
CN| : compositional variety. In this paper we develop models for fluids where the variety of constituents arises from an 
^ internal symmetry group. We envision applying our theory to high densities of non-Abelian quarks, with or without 

, additional interaction to a dynamical gauge field. 
^ ^ ■ High energy collisions of heavy nuclei can produce of a plasma state of quarks and gluons. This new state of matter 
' has recently been of great interest both theoretically and in experiments at the RHIC facility and at CERN. Most of 
I the theoretical work in this area has been based on perturbative Quantum Chromodynamics at high temperatures with 
hard thermal loop resummations and other improvements Q. This can be a good description at high temperatures 
and for plasma states that are not too far from equilibrium. An alternative approach, which may be more suitable 
' for nondilute plasmas or for situations far from equilibrium, would be to use a fluid mechanical description. 
, It is well known that the equations of fluid mechanics can be derived from particle dynamics by taking suitable 
• averages of Boltzmann-type equations. Specifically for the quark-gluon plasma, some work along these lines was done 
' many years ago using single particle kinetic equations Q. These kinetic equations form a hierarchy, the so-called 
BBGKY hierarchy, involving higher and higher correlated iV-particle distribution functions. To be able to solve 
' them, one needs to truncate the hierarchy, very often at just the single particle distribution function. Therefore, 
the feasibility of solving these equations limits the kinetic approach to dilute systems near equilibrium, where the 
truncation can be justified. We might therefore expect that the regime of validity of equations derived within kinetic 
theory is likewise limited. However, fluid dynamical equations can also be derived from very general principles, 
showing that they have a much wider regime of validity, and, indeed in practice, we apply them over such a wider 
range. This is the 'universality' of fluid dynamics. 

In the context of a non-Abelian plasma, in analogy with the ordinary fluid mechanics, we may therefore ask for an a 
priori derivation of a non-Abelian fluid mechanics, which incorporates the non-Abelian degrees of freedom, coupling to 
a non-Abelian gauge field, etc. This theory may be valid for dense, nonperturbative and nondilute systems. Further, 
a canonical or symplectic formulation (at least in the conservative limit) is important for quantization. At the same 
time, the analysis based on the kinetic equations remains useful to us as a guide for arriving at Lagrangian, canonical 
description. 

Our goal is to provide plausible equations that generalize to the non-Abelian situation the continuity 

9tp(t,r)-f V- (p(i,r)v(t,r)) =0 (1) 

and Euler equations 

(9tv(t, r) + V • Vv(i, r) = force (2) 

which govern conventional Eulerian perfect fluids. In (^ force denotes forces acting on a unit volume of the fluid. 
The above four equations or their relativistic versions can be equivalently presented as the four conservation laws for 
the energy and momentum densities (either nonrelativistic or relativistic). 
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In the non-Abelian generalizations one must first fix the group transformation law for the collective variables and 
the covariance properties of the equations, which now may number more than four, and cannot be comprised solely 
in energy-momentum conservation. An obvious additional equation is the (covariant) conservation equation for the 
non-Abelian current, but the relation of this current to the collective variables needs to be established. 

Our constructions begin by postulating collective variable Lagrangians that are invariant against space-time trans- 
lations and internal symmetry transformations. Thus energy-momentum conservation and current conservation are 
assured by Noether's theorem, while the Euler-Lagrange equations describe the temporal evolution of the collective 
variables. In the Section |l| we review the Lagrangian as well as other aspects of nonrelativistic and relativistic Abelian 
fluids. This is done to motivate and contrast our non-Abelian constructions, which we develop in Sections lU and IV. 
Derivations are postponed to the Appendices; in the text only results are presented and discussed. 



II. REVIEW OF ABELIAN FLUIDS 

Equations and (|^) can be found as the Euler-Lagrange equations for the Lagrange density 

1 
2 

where our coordinates are x'^ — {ct, r). The various quantities are defined as 



C = -ra^ + -p^r^ - V (3) 



f = (cp, pv) (4) 
flp = 9^6' + adfj^P (5) 

and is a p-dependent potential that gives rise to the force. (In the usual theory, V depends only on p, but we 
allow dependence on derivatives of p as well.) The canonical 1-form is determined by the first contribution to C while 
the dynamics is encoded in the ^pv^ — V term. Specifically, variation of v shows that v is given by the Clebsch 
parametrization jj] 

V = + aV/? . (6) 

Variation of results in (Q), while variation of the Clebsch potentials (a, (3) produces subsidiary equations 

j/^a^a = p [dta + V • Va) = , . 

fd^l3 = p{dt(3 + v-Wf3)^0 ^'^ 

which are needed in the subsequent derivation of the Euler equation. Finally, the "Bernoulli" equation emerges by 
varying p 

dt9 + adt(3+^^-^J drV . (8) 

( J- / drV is just the Euler-Lagrange derivative of V .) The Euler equation (||) then follows by taking the gradient of 
and using (j^). The force on the right of (j^) is now seen to be 

force = -V^ [ drV (9) 
5p J 

The energy £ and momentum V densities carried by the fields are 

£ = ipv2 + y (10) 

V = pw . (11) 

These satisfy continuity equations as a consequence of (or equivalent to) (|l|) and (||). 

Interaction with electromagnetic fields is accommodated by including the Maxwell Lagrange density and adding 
the electromagnetic vector potential to a^, whereupon becomes the electromagnetic current. This gives rise to the 
Lorentz force in the Euler equations, thereby producing magnetohydrodynamics. Note that in this model the mass 
density moves with the same velocity as the charge density as is seen by inspecting V and j. 
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An interesting realization of the above equations is provided by Madelung's "hydrodynamical" rewriting of the 
Schrodinger equation Q 



ifidtip{t,v) 



2m 



where the wave function is presented as 



(12) 



(13) 



The imaginary part of (|T^) reproduces the continuity equation (^, when V9 is identified as the velocity v, with 
vanishing vorticity V x v = 0. Also the quantum current ^liiiip*Vip become pv. The real part of (|l^) gives the 
Bernoulli equation with a and /3 set to zero and with a "quantum" force derived from 



The Euler equation (||) follows by taking the gradient of the Bernoulli equation. 
The Lagrange density for a relativistic fluid is chosen as |3| 



(14) 



^ = -i^flp - /(") ' ^ = \Af%A? ■ 

The canonical 1-form is as in the nonrelativistic case (|)-(|) and now is written in the Eckart form 

u^ = (c,v)/v/r^W?, 



p\/ 1 — jc?- 



(15) 



(16) 



/(ri) encodes the dynamics; for the free theory f(n) = nc? . Variations of 0, a and /3 still produce equations (|i|) and 
(0) while variation of evaluates 



r.2 



fin) 



(17) 



where dash denotes differentiation with respect to the argument of the function. We call (^^ the relativistic Bernoulli 
equation, which leads to the Euler equation with the following steps. The curl of (|l^) reads 

(18) 



dt,a„ - d^a^ = \ -^f {n)j ~ y-^f (n) 

Since according to (^ the left side is equal to d^ad^P — d^ad^,[}, it vanishes by (Q) when projected on u^. Thus there 
remains the relativistic Euler equation 



It is easy to show that in the nonrelativistic limit, with 



n 



(c,v) 



and 



f = nc^ + V{n) 



the Lagrange density (|l|) passes to (||) (apart from an irrelevant term —pc 
reproduces the Euler equation (^. 

The energy-momentum tensor for ( |l5| ) reads, after is eliminated with (17) 



(19) 

(20) 
(21) 

(22) 

and the spatial component of ( |l9| ) 



TA'^'^-^M^ („/'(„)_/(„)) + 



■nfin) 



Its divergence d^T'^'^ can be expressed as 



— a^(«V'("))-9V'(^) 



(23) 



(24) 
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The first term on the right vanishes by the virtue of the continuityequation (|l|) and the second vanishes by the 
relativistic Euler equation ( |l9| ) ; or vice versa: since the two terms in (|2^ ) are linearly independent (the first is parallel 
to u'^ and the second is orthogonal to it), conservation of T^" implies (|l|) and (|l9|). Relativistic magnetohydrodynamics 
is achieved, as in the nonrelativistic case, by adding the electrodynamical potential to a'^. 

Finally we record one fact, which we shall use below, about the alternative, Lagrange, formulation of fluid mechanics. 
Here one describes the fluid with co-moving coordinates X(t, x) where x is a (continuously) varying particle label. It 
may be chosen to coincide with X at i = 0. The relation between Euler and Lagrange variables is the following. For 
the density we have 

p(t, r) = J dxS (X(t, x) - r) (25) 

which is normalized to unity at i — 0. (The integral and the (5-function follow the dimensionality of space.) Evidently, 
in the course of the x integral, the J-function evaluates x at an expression x(^i^) such that X(i, x(t,r)) = r; i.e., x 
is the inverse of X. There is also a Jacobian. Thus 



1 , dX 
- = det 



(26) 



The Euler velocity is given by 



or 



p dx^ 

v=5,XU^ (27) 

p{t,r)v{t,r) = J dxdtX{t,yi)S{X{t,x)-r) . (28) 

A simple calculation shows that the density p and current j = pv defined by these equations satisfy the continuity 
equation (|l|). [The Euler equation follows by differentiating ( |28| ) with respect to time, positing a force that determines 
d^X/dt^, and performing the corresponding x integral; we shall not need this here.] 



III. NON-ABELIAN MODELS BASED ON A PARTICLE SUBSTRATUM 

A. Non-Abelian current 

Before presenting a specific model, we give a general analysis of the non-Abelian current J^' — (cpa,Ja)- 
The conventional formula for the current of a single, non-Abelian point particle, moving in a 4-dimensional space- 
time {t,r}, along a space-time path X'^{t) (t parametrizes the path) is 

dr Q,{r)^^5{X\r) - ct)S{X{r) - r) . (29) 

This is covariantly conserved 



proved that Qa satisfies the Wong equation 

dQair) 



d^^JH + fabcAU^ = {D^J^^)a = (30) 



+ /a6c-^X^(T)A':(X(r))Q,(r) = (31) 



For several particles Qa and X^ are indexed by a discrete particle label n, which is summed in the definition of the 
current Jj^. In a continuum limit n — > x we have 

dx dr Qair, x) ^ ' s (X"(t, x) - ct) S (X(t, x) - r) (32) 

^Qa(r,x)-|-/,b,— Xnr,x)4(X(T,x))Q,(r,x) =0 (33) 
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[Notice that, since we replace a sum over n with an integral over x, Qa{T, x) is now a charge density.] The parametriza- 
tion may be fixed at X°(t, x) = cr so that the equations read 



Pa{t,r) = J dxQa(i,x)5(X(i,x)-r) 
3a{t,r) = [ dxQait,^)dtX{t,K)6{X{t,K)-r) 



(34) 



dtQait, x) + fabc [cAUt, X(t, x)) + dtX{t, x) • A''(^, X{t, x))] Q,{t, x) = 



(35) 



Observe that just as in the Abelian case discussed above, the x integration evaluates x at x{t,J^), the inverse of X, 
and the Jacobian factor (^detdX^dx^) ^ is just the Abehan charge density p [see ( p^ ) ]. Thus 



Pait,r) = Qa{t,r)p{t,r) 
Ja(i,r) = Qa{t,r)p{t,r)v{t,r) 



(36) 



j;^{t,r)^Qa{t,r)r{t,r) , (37) 

where 

Qa{t,r) = Qa(i,x)|,^^ (38) 

p{t,r)Qa{t,r) = y da:Q,(t,x)5(X(t,x)-r) (39) 

Moreover, differentiating ( ^9|) with respect to time and using (|l|) and ( |35| ) results in an equation for dQ'^/dt, 

dtQait, r) + v(t, r) • VQait, r) - -/,bc (cA|;(t, r) + v(t, r) • A'it, r)) Q,(t, r) (40) 
which can also be written as 

r{D^Q)a = . (41) 

This is analogous to the Abelian equations (0). Equations (^0|) and (|4^) can be understood from the fact that the 
covariantly conserved current ( |30|) factorizes according to (|37^ into a group variable Qa and a conserved Abelian 
current j^^. Consistency of (|^), (|30| ) and (37|) then enforces (|4l|). 

We recognize that formulas (p^ and (3#) are the non- Abelian version of the Lagrange variable - Euler variable 
correspondence [see ( p5| ) and (|28D]. Also, (3^, ( |4C| ) and ( |4l| ) are the field generalizations of the particle Wong equation 
(|3T|), with ( |35| ) being presented in the Lagrange formalism and (40), ( |4l| ) in the Euler formalism. The decomposition 
of the non- Abelian current in ( |37| ) is the non- Abelian version of the Eckart decomposition ( p^ ) [^. Indeed, ( p7| ) may 
be further factored as in ( |l^ ) 



J,^(t,r) = Q,(t,r)n(t,r)ii^(t,r) , 



(42) 



In the remainder of Section III, we are guided in our construction of a dynamical model for non- Abelian fluid 
mechanics and "color" hydrodynamics by the above properties of the non- Abelian current, which follow from the very 
general arguments, based on a particle picture for the substratum of a fluid. In Section ^ we present a different 
model, based on a field theoretic fiuid substratum. 



B. A Model for non- Abelian Color Hydrodynamics 

The model that we present is based on a group with group elements g, and anti-Hermitian Lie algebra elements 
with basis T". 



j^/j^a jabcrpc 

tr {T'^T'') = -^S"'' . 



(43) 
(44) 
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The Lagrange density for an Eulerian fluid built on such a group is taken to be the following generalization of the 
Abelian expression ( |l5| ) 

C = f2 tr [Kg-^D^g] ~ f{n) + Cgauge ■ (45) 
Here is an Abelian vector field (current) which can also be decomposed as in (|6|) 

= {cp, pv) — nu^ , u'^u^ = (46) 

The covariant derivative 

D^^g - d^g + A^g (47) 

involves a dynamical non- Abelian gauge potential yl^ = v4°Ta whose dynamics is provided by Cgauge- K is a fixed, 
constant element of the Lie algebra. The first term in C contains the canonical 1-form for our theory and determines 
the canonical brackets, while /(n) describes the fluid dynamics. The theory is invariant under gauge transformations 
with group element U 



9 U-^g 



-If A , f^^rr (48) 



According to (45), the current J^' to which A'^ couples is of the Eckart form ( [4 2D 

= (49) 

with 

Q = QaT'' = gKg-^ . (50) 

For consistency of the g auge field dynamics J'^ = J^T" must be covariantly conserved. To satisfy the general 
arguments of Section III A, j'' must be divergenceless and then the Wong equation ( ^l| ) will follow. In Appendix ^ we 
show that both conservation laws are a consequence of invariance of the action with respect to variations of the group 
element g; arbitrary variations of g lead to covariant conservation of J'' (^0|) while the particular variation Sg = gK\ 
ensures that is conserved (||). Therefore the Wong equation (^ is also a consequence. 

Thus our model reproduces all the equations satisfied by the current that were established in Section [II A from 
general principles. Indeed the canonical (first) term of the Lagrangian (^) is like a Kirillov-Kostant 1-form, which 
has been previously used to give a Lagrangian for the point particle Wong's equation (^ij) |9|. Moreover, as we show 
in the Appendix H, the canonical brackets implied by the canonical 1-form ensure that the charge density algebra is 
represented canonically 

{pa{t, r), pb{t, r')} = fabcPcit. r)5(r - r') . (51) 
It remains to derive the Euler equation. This is accomplished by varying j'' ; stationary variation requires 

2ir [Q{D^g)g-^]^^r{n) (52) 



which we call the non- Abelian Bernoulli equation. The Euler equation then follows, as in the Abelian case, by taking 
the curl 

(2 tr [Q{D,g)g-'] ) - 9. (2 tr [Q{D^g)g-^] ) = (^/'(")) - (^/'(^)) • (53) 
In the Appendix [c| we show that manipulating the left side allows rewriting (53) as 

2 tr [{D^Q){D,g)g-^] + 2 tr [QF^,] = 8^ (^/'(")) " (^/'H) ■ (54) 
Finally, contracting with — nu^ and using ( ^l| ) produces the relativistic, non- Abelian Euler equation 

—d^ {u,,f'{n)) - nd,f'{n) - 2 tr [J'^F^,] . (55) 
The left side is as in (^9|) while the right side describes the non-Abelian Lorentz force acting on the charged fluid. 
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Apart from the gauge field contribution, tlie energy-momentum tensor for the Lagrange density ( ^ ) is the same as 
in ( p3| ) when ( p2| ) is used to eliminate tr Kg~^D^g 



T^- = -g^^ (nfin) ~ fin)) + ^nf'{n) 



(56) 



Just as in the Abelian case, the divergence of T'^'' entails two independent parts: one proportional to u'^ and the 
other orthogonal to it 



W-f'in) 



The first vanishes by the virtue of (Jl]) and the rest is evaluated from Euler equation (|55|), leaving 

d^T^''' = 2 tr [J^F^""] 
which is canceled by the divergence of the gauge-field energy-momentum tensor 

'^tJ^'^gauge — ~2 tr [J fi,F^ ] . 



(57) 



(58) 



(59) 



Thus energy-momentum conservation reproduces Abelian current conservation 1l) and the non-Abelian Euler force 



equation (^6|), but the Wong equation ( [flD h as to be enforced additionally. This is achieved by our Lagrangian (|45|). 

We record the nonrelativistic limit of (|55D; using ( pO| ) - (p2[), we find that the nonrelativistic limit for the spatial 
component of (|5^) gives the Euler equation with a non-Abelian Lorentz force 



5tv + V • Vv = force + Q"E„ 



(60) 



where force is the pressure force coming from the potential V (and is therefore Abelian in nature), while non-Abelian 
force terms involve the chromoelcctric and chromomagnetic fields 



c 

_ pa 



(61) 



It is seen that our non-Abelian fluid moves effectively in a single direction specified by j = pv. Nevertheless, it 
experiences a non-Abelian Lorentz force. In the next subsection we present a generalization wherein the non-Abelian 
fluid develops several independent directions of motion. 



Generalization 



A generalization of our Lagrange density (45) that will give rise to several fluid components carrying various 
densities and moving with various velocities is obtained by choosing several directions in the Lie algebra, -ft'(s), and 
coupling to different Abelian currents M 



^ = ^i(s)^^'^ [K(s)9 ^D^,g]- f{n(i),n,^2),---,n(^r))+^gauge 

Jis) = (cP(s),P(.)V(,)) (62) 

(Sums over s are indicated explicitly; the summation convention does not apply.) Evidently the current which couples 
to the gauge potential is now 

r 

J'^^Y. Q(s)^(s) ' with Q(,) = A'(,)5-i . (63) 

s=l 



Arbitrary variation of g ensures that (|63| ) is covariantly conserved, but we also need the conservation of individual 
j^jj. This is achieved by choosing special forms for 5[s)g = gK(^s)^{s) which will work provided that different K(^g-j 
commute (see Appendix Therefore, we choose the K(^g^ to belong to the Cartan subalgebra of the Lie algebra and 
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the total number of different channels equals the rank r of the group. With this choice for the i^(s), special variations 
of q ensure 



^mjII.) = dtp(^s) + V • (p(s)V(s)) = 
The Wong equation which follows from the conservation of the non-Abelian current now reads 

r 

E^T.)^mQ(.)=0. 

5-1 

Varying the individual produces the Bernoulli equations 



2tr [Q(,)p^5)g-i] = ^/(^) where /< 



9n, 



is) 



-/("-(l),"-(2), ■ • • ,"-(r)) 



Again, the curl of the above can be cast in the form 

^ H)/^^') - } = 2 tr [P^Q(,))(i?''5)g-i] + 2 tr [Q(.)f^^1 



When contracted with j 



is) = 



this leaves 



(64) 



(65) 



(66) 



(67) 



(68) 



However, unlike in the single channel case, the right side does not simplify since j^(s)Q{s) cannot be replaced by 
because the latter requires summing over s. Also the first right-hand term in (^) does not vanish since ( |65| ) requires 
summation over s. 

The energy-momentum tensor is 



Its divergence of course reproduces (57) 



+ n 



(69) 



(70) 



The first term in the brackets vanishes according to (M) and the remainder is evaluated from (pq) as 



r 



Since now we are summing over all channels, it follows from (|63| ) and (65) that, as before, 

d^.T"" 2 tr [ J^i^'"'] 

A more transparent picture of what is happening is given if the dynamical potential separates 

r 

/(n(i), . . . ,n(^)) ^ f{s){n(s)) 



1 

is) 



(71) 

(72) 
(73) 



Then the left side of (|68| ) refers only to variables labeled s, while the right side may be rewritten with the help of (|6: 
and KB) to give 



r 

d, - ^is)d''fls) = 2 tr [J,Fn - 2 E tr [j^^,,) {Q(s-)F^'' + {D^^Q(s-)){D'' 9)9-')] ■ (74) 
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Thus in the addition to the Lorentz force, there are forces arising from the other channels s' 7^ s. Note that for 
separated dynamics ([72|), the energy-momentum tensor also separates 



= = y: -3- (-(.)/('.) - /(.)(-(.))) + ^^-(.)/(.) 



(75) 



but the divergence of individual T^'^^ does not vanish. This reflects that energy is exchanged between the different 
channels and with the gauge field, as is also evident from the equation of motion (^^. ft is clear that this fluid moves 
with r different velocities V(s). 

The single-channel Euler equation (55) is expressed in terms of physically relevant quantities (currents, chromo- 
magnetic fields); the many-channel equation (68) involves, additionally, the the gauge group element g. One may 
simplify that equation by going to special gauge, for example g — /, so that the right side of (pq) reduces to 



2tr [j^is)iD^^Qis)){D,g)g-']+2tv [j(%Q(,)F^. 
while the Wong equation (|65|) becomes 



= 2tr 



(76) 



{s)\ 



(77) 



It is interesting that in this gauge the non- linear terms in F^'^ disappear. 



IV. NON-ABELIAN FLUIDS WITH A FIELD SUBSTRATUM 

In the Abelian case, the Madelung parametrization ([l3| ) of the Schrodinger equation gives the conventional 
nonrelativistic Euler equation, [even while the forces are derived from a potential that depends on density and 
(unconventionally) on its derivatives] (see Section ||) . We are therefore led to examine a Madelung-like construction 
for a non- Abelian, "colored" Schrodinger equation 

indt^p = -—V^i: (78) 
Zm 

We consider only the free, nonrelativistic case, and the "non-Abelian" structure resides solely in the fact that the tp 
is a multi-component object, transforming under some representation of a group. The color degrees of freedom also 
lead to the conserved color current 

J^^{cpa,3a), Pa=i^^T''iP, J„. = -Re V^T" • (79) 

rn 



Of course, the singlet current 



r^icp.j), P = ^H: j--Im^tv^. (80) 

ni 



is also conserved. For definiteness and simplicity, we shall henceforth assume that the group is SU{2) and that the 
representation is the fundamental one: = a°'/{2i), {T°',T^} = —6°'''/2. We shall also set the mass m and Planck's 
constant h to unity. The non-Abelian analog of the Madelung decomposition (|l^) is 

= VP9u (81) 
where p is the scalar V'^V': 5 is a group element, and m is a constant vector that points in a fixed direction [e.g., for 
Q ^ , then iu^T°-u ~ 5°''^ /2\. The singlet density is p, while the singlet current j is 

j = pv , V = -iu^ g-^S/gu . (82) 
With the decomposition (|8|), the color density ( |7^ ) becomes 

Pa^QaP, Qa^iu^ g-^T''gu^iRabU^T^u= Rabt'' /2 (83) 
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where Rat is in the adjoint representation of the group and the unit vector t° is defined as t°'/2 = iu^T°'u. On the 
other hand, the color current reads 



which with the introduction of 



becomes 



V = v'^t" 



P 



(84) 

(85) 
(86) 

(87) 



Unhke the Abehan model, the vorticity is nonvanishing 

V X v'' = e"''^v^ X 

A difference between the Madelung approach and the previous particle based one is that the color current is not 
proportional to the singlet current. Equation (p7\) may be written as 



where the "orthogonal" velocity v° is defined as 



4 = (,5'''' - ft'') v' 



(89) 



(90) 



Equation ( p8| ) shows that color current possesses components that are orthogonal to the singlet current. 

In Appendix |^ we derive for the SU (2) case the decomposition of the Schrodinger equation with the parametrization 
(|8T|). Two equations emerge: one regains the conservation of the Abelian current (|l|) and the other is the "Bernoulli" 
equation 



{9~'d,gr 



VP 



e + -V • (pe''''Vt'=^ 

p 



(91) 



It is further verified that the covariant conservation of the color current is enforced by both (|l|) and (^l|). However, 
there is no Wong equation because the color current is not proportional to the conserved singlet current. Finally, 
using the identity, which is a consequence of the definition (851) 



:V 3-194.9 +6""^" (g-^dtg 



one can deduce an Euler equation for 9tv° from (^. 
We record the energy and momentum density 



(92) 

(93) 
(94) 



Both parallel and orthogonal components of the velocity contribute to the energy density but only the parallel 
component v contributes to the momentum density. It is clear that within the present approach the fluid color flows 
in every direction in the group space, but the mass density is carried by t he un ique velocity v. This is in contrast 
to our previous approach where all motion i s in a single direction (Section III B ) or at most in the directions of the 
Cartan elements of the Lie algebra (Section tllC). 

The difference between the two approaches is best seen from a comparison of Lagrangians. For the color Schrodinger 
theory in the Madelung representation 



i 1 



ipu^ g ^dtgu- -pv° • v'' 



(Vpl 
8p 



(95) 
(96) 
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With u ® =1 12 — 2iK, the free part of the above reads 



J^lchrod^nger = P 2 tr [Kg~'dtg] - -pv-^ • v'^ (97) 



On the other hand, the free part of the Lagrange density (Hsh in the nonrelativistic hmit is 



£" = p2tr [Kg-^dtg] + pv • 2 tr [Kg^'^Vg] - y^p^ic^ -v'^) 
« p 2 tr [Kg-^dtg] + pv • 2 tr [Kg-^Wg] - pc^ + ^pv^ 

= p 2 tr [Kg-'dtg] - \pv^ - pc^ (98) 

where we have used v = —2 tr [ivr(7~^V(7] , which foUows upon the variation of v, in the next-to-last equahty above. 
Thus the canonical 1-form is the same for both models while the difference resides in the velocity dependence of their 
respective Hamiltonians. Only the singlet v enters (^) while the Madelung construction uses the group vector v"^. 

Finally, note that while the Euler equation, which emerges when ( ^l|) and (^2|) are combined, intricately couples all 
directions of the fluid velocity v*^, it does admit the simple solution v° = vt°, with v obeying the Abelian equations 
that arise from (|l2|)-(p|). 



V. DISCUSSION 



In this paper, we have presented in Section III two distinct non-Abelian generalizations of ordinary particle based 
Abelian fluid dynamics. Both versions use a fluid generalization of the Kirillov-Kostant form that naturally encodes 
the algebra of charge densities dsi]), which is need ed if the charge density is to be identified as the gene rator of 



non-Abelian symmetry transformations. In Section [II C we generalized the first version, given in Section IIIB so 



that the density is specified in terms of a set of Abelian densities equal in number to the rank r of the Lie algebra, 
rather than a single density as in the first version. Since the charge density at a point in the fluid is an element of 
the Lie algebra, diagonalization shows that an invariant specification must use r eigenvalues. Alternatively, we may 
use the r Casimir invariants of pa at a point to characterize it. Therefore the appearance of r Abelian currents in 
the Lagrangian is entirely natural. This is also in accord with what happens with the Kirillov-Kostant form where 
our j'l'^j are replaced by the fixed weights of a representation of the Lie algebra and lead to that representation upon 
quantization. 

The two versions also differ in the formula for the current. Our first version, which is mathematically more concise, 
gives the non-Abelian Eckart decomposition, (^9|), while the second version does not allow the factorization of the 
current into a non-Abelian charge density and an Abelian velocity, rather it is a sum of such factorized expressions - 
a generalized Eckart decomposition as in ( |50| ). The Eckart decomposition shows that we can choose a local Lorentz 
frame for which a given fiuid element can be brought to rest. The charge density is then related to the charge carried 
by this element. By contrast, in the second version, with a generalized Eckart decomposition, we see that even if we 
choose a frame where one Lie algebra component of the velocity is zero, the other color velocities need not be. Thus 
the latter applies to a situation where color separating flows can occur. Physically, it is not yet clear what kinematic 
regimes of a quark-gluon plasma, for example, would admit or require such flows. 

We note that the currents we have obtained are the non-Abelian analogues of the irrotational part of the Abelian 
current, even though the vorticity is nonvanishing. The other components can be easily incorporated, if needed , by 
generalizing the Lagrangian in (|62|) as 



= ^ {2tr [K(^s)9 ^£^^5] +aA<(s)} -/("(i),"(2),---,".(r)) + ^SQ«se (99) 

s=l 



where a^(s) is given by 



ap(s) = (100) 

The final fluid equations remain unchanged. 



We have also derived in Section IV a field based fiuid mechanics by extending the Madelung construction to the 
non-Abelian situation. Here the "Euler" equations are much less appealing because they involve velocities in all group 
directions. We know of no compelling physical reason for preferring this field based model over the particle based 
one, even though in the Abelian case it coincides with the particle-based construction. 
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APPENDIX A: VARIATIONS OF g 

We determine the variation of 

/O = / dtdrJ2 J is) 2 tr K^,^g-^D^g (Al) 

when g is varied either arbitrarily or in the specific manner 

g-'Sg^K^,,~,X (A2) 

where A is an arbitrary function on space-time. This will provide the needed results (|[) and ( |30| ) for the single channel 
situation as well as (^J) and (^5|) for many channels. 

Recall the definitions Q(^s) = 9^(3)9^^ ^-nd Dfj^g = d^g + A^g, which implies D^g^^ = df^g^^ ^ ff^^^P- First, the 
variation of g^^Z^^g is established 

S {g-^D,,g) = ~g-Hgg-^D^g + g-^D,,5g (A3) 
To evaluate the last term, note that D^5g — {gg^^Sg^ — {Df^g)g^^5g + gD^ (^g^^Sg). Thus 

5 (g-'D^g) = (g-'Sg) + [g'^D^g, g-'6g] . (A4) 
Inserting (A4) into the variation of 7° in ( |AlD , integrating by parts, and rearranging the trace with ivT^^) gives 



/r 
dtdrJ2{d^jis)^^^Kis)g-'Sg + j^^^2tT [g-'D^g,Kis)]g-'Sg) . (A5) 
s=l 

Considering first arbitrary variations: the vanishing of (5/*^ requires 

r 

E {^Jis)Kis) + Jis) [9-'D,g, ) = (A6) 
or, after sandwiching the above between g . . . g~^, 

r 

E + [{D,g)g-\ Qis)]) - . (A7) 

s=l 

Finally we verify that 

[^p.93"\0(.)] , (A8) 



SO that the desired resuhs (30) and (|65|) foUow 



r / \ 



Next we consider the specific variation (A2) and separate the sum (KS) into the term s — s' and s ^ s'. After a 



rearrangement of the last term in ( A£ ) , we get 

«° = - y dtdr(a^j^^,^2tri^(,,)if(,,)A + j^^,^2tr.9-ii?^g[i^(,,),if(,,)]A) 

+ E 2 + ^r^) 2 tr g-'D^g[Kis),Kis')]x) (AlO) 

The first commutator vanishes; so does the second when K(^s) and K(s') commute, i.e., when they belong to the 
Cartan subalgebra. Also 2 tr Kf^g^K^^gi^ — --K^^^K^^,y, for s' — s this is constant, while for s' 7^ s it vanishes when 

it is arranged that distinct elements of the Cartan algebra are selected. Thus for stationary variations j^^.^ must be 

conserved, and (III) as well as (|30|) are validated. 
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APPENDIX B: CHARGE DENSITY ALGEBRA 



The portion of the Lagrange density ( p5[ ) that determines the Poisson bracket is 

Ccanomcai = p2tr Kg^^dtg = p2 tr Qdtgg^^ ■ 



(Bl) 



With a parametrization of the group element, e.g., g{(p) = "^"j one sees that (dtg) g ^ has the form —9(iy9aC'^((,c)r'', 
where the non-singular matrix C1{ip) is defined by 



Thus 



and the momentum conjugate to (pa is 



O^Pa 

^canonical — pdtPaC%Qb = dtipaClpb 

= cipb . 



(B2) 

(B3) 
(B4) 



Pa = c^^n" 



(B5) 



The non-Abelian charge density pa is a function of (t, r) and for (51) we need the bracket with another density 
evaluated at {t,r'). Since the dependence of c^j, on ip involves no spatial derivatives of p, it is clear that the brackets 
will be local in r — r'; just as is the bracket between p and 11. 



{pa{t,rhpb{t,r')} = (c\,^U'^' - a^b] 5 {r-r') 



-cVct.,^^c'=„,n" - a ^ H 5 (r - r') 



dpb' 



-b"-c' 



dpb' 



a ^ 6 (5 (r — r') 



(B6) 



To evaluate the derivative with respect to p, return to (B2) and observe 



9 ^2trS^g-^r" 



dipb' dpb' \ dp 
9^9 



2tr 



dpc'dipc 



(B7) 



The first term in the parentheses is symmetric in (&', c'); when inserted in (B6) it produces a symmetric contribution 
in (a, 6) and does not contribute when antisymmetrization in (a, 6) is effected. What is left establishes (|45|). 

{pa{t,v),pb{t,v')} = (c\,c^,,C'^'^,C''^.2irT'''T''''r'' p,n - a ^ b) 5 {v ~ v') 

= -{2iY T^'T^T''" p^n - a^b^5{v-r') 

= -2ti r'^^T'^T"" p^„5{r-v') 

= r"'Pc{t,v)6{v-v') (B8) 



APPENDIX C: MANIPULATING EQUATION (|53D 

Observe that the first term in ( |53| ) equals 

9^2 tr [Q{D,g)g-^] = 2 tr [{DMDu9)9~^ + Q{D^D,g)g-^ ~ Q{D,g)g-\D^g)g-^ .] 



(CI) 
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The first term on the right side is rewritten with the help of ( [A^ ) and combined with the last term, leaving 

2tr {Q{D^D,g)g-^ ~ Q[D^g)g-\D,g)g-^) 
After antisymmetrization in (/i, J^), the left side of ( p3| ) reads 

2 tr Q {{[D^, D,]g)g-^ - [{D^,g)g-\ {D,g)g-^]) - 2 tr [QF^, - [Q, {D ^g)g-^]{D,g)g-^) (C2) 



When (A8) is used again, (pa) becomes the left side of (p 



APPENDIX D: NON-ABELIAN MADELUNG PARAMETRIZATION 



When ( ^iD is inserted into ([78|), and use is made of the definition (p5[), we find in the SU(2) case 

^zdtpu + ipig-'dtgrr- u = ~\^py^^pu + V(pv'^)T'^ u + ipv-^ • v"^ u (Dl) 



Next (Dl) is prcmultiplied by u^, where it implies 

idtp + p{g-^dtgTe = - VPV^ VP - *V (pv'^t'^) + pv"^ • v"^ . (D2) 



The imaginary part reproduces the continuity equation for the singlet current (^2|) , while the real part gives 

{g-^dtgYf^ = -^^^^p + ^^ -w^ . (D3) 



To obtain further information, we premultiply ( |DlD with u^T^ . This gives 

dtpt" - ip{g-'dtgf + pig-'dtg)\'"'H- = i^p^^^pt" ~ V • (pv") - zV • (pv"^) e'^-t^ - *pv'^ • v'^t" . (D4) 
The imaginary part gives (|9^) while the real part is identically satisfied by virtue of (|l|) and (^l|). 
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